Phonon-induced decay of the electron spin in quantum dots 
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We study spin relaxation and decoherence in a GaAs quantum dot due to spin-orbit interaction. 
We derive an eflective Hamiltonian which couples the electron spin to phonons or any other fluc- 
tuation of the dot potential. We show that the spin decoherence time T2 is as large as the spin 
relaxation time Ti, under realistic conditions. For the Dresselhaus and Rashba spin-orbit couplings, 
we find that, in leading order, the effective magnetic field can have only fiuctuations transverse to 
the applied magnetic field. As a result, T2 = 2Ti for arbitrarily large Zeeman splittings, in contrast 
to the naively expected case T2 <C Ti. We show that the spin decay is drastically suppressed for 
certain magnetic field directions and values of the Rashba coupling constant. Finally, for the spin 
coupling to acoustic phonons, we show that T2 = 2Ti for all spin-orbit mechanisms in leading order 
in the electron-phonon interaction. 
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Phase coherence of spin in quantum dots (QDs) is of 
central importance for spin-based quantum computation 
in the solid state P, 01 • Sufficiently long coherence times 
are needed for implementing quantum algorithms and 
error correction schemes. If the qubit is operated as a 
classical bit, its decay time is given by the spin relax- 
ation time Ti, which is the time of a spin- flip process. 
For quantum computation, however, the spin decoher- 
ence time T2 — the lifetime of a coherent superposition of 
spin up and spin down states — must be sufficiently long. 
In semiconductor QDs, the spin coherence is limited by 
the dot intrinsic degrees of freedom, such as phonons, 
spins of nuclei, excitations on the Fermi surface (e.g. in 
metallic gates), fluctuating impurity states nearby the 
dot, electromagnetic fields, etc. It is well known (and 
experimentally verified) that the Ti time of spin in QDs 
is extremely long, extending up to 100 /is. The decoher- 
ence time T2, in its turn, is limited by both spin-flip and 
dephasing processes, and can be much smaller than Ti, 
although its upper bound is T2 < 2Ti . Knowledge of the 
mechanisms of spin relaxation and decoherence in QDs 
can allow one to find regimes with the least spin decay. 

Recently, the spin relaxation time Ti in a single- 
electron GaAs QD was measured by means of a 
pulsed relaxation measurement technique (PRMT) 0. 
This technique was previously applied to measure the 
triplet-to-singlet spin relaxation in a two-electron quan- 
tum dot 0, yielding a spin relaxation time of 200 /xs. 
The application of the PRMT to Zeeman sublevels be- 
came possible with resolving the Zeeman splitting in dc 
transport spectroscopy 0, U , which required a magnetic 
field S > 5 T. The results of Ref. show that Ti > 50 ^is 
at B = 7.5 T and 14 T, with no indication of a i?-field de- 
pendence. Experimental values for the spin T2 time in a 
single QD are not available yet, but an ESR scheme for its 
measurement has been proposed The ensemble spin 
decoherence time T2 was measured in n-doped GaAs bulk 
semiconductors 0, demonstrating coherent spin preces- 



sion over times exceeding Tj* ^ 100 ns. This indicates 
that the decoherence time of a single spin is even larger, 
T2 > Tj*- However, the mechanisms of spin decoherence 
for extended and localized electrons are rather different 
(cf. Ref. •?) and, therefore, the result of Ref. Igl is not 
conclusive for electrons in QDs. 

Different mechanisms of spin relaxation in QDs have 
been considered, such as spin-phonon coupling via spin- 
orbit 0, ^3 oi' hyperfine interaction jljj, and spin- 
nuclear coupling [13, uM, ■ The spin-orbit mechanisms 
yield no spin decay at J5 = 0, due to the Kramers degen- 
eracy of the dot states. Interestingly, for GaAs QDs, the 
orbital effect of B leads to n o sp in decay in lowest order 
in spin-orbit interaction |lOl llll II6L [iTf . This is due to 
the special form (linear in momentum) of the spin-orbit 
coupling in two-dimensions (2D). The leading order con- 
tribution is, thus, proportional to the Zeeman splitting 
and leads to long spin-flip times Ti in GaAs QDs varying 
strongly with the B-field j^. However, previous theo- 
ries do not apply to the high values of B used in recent 
experiments pli ^'^d thus, no comparison could be made 
so far. As for the nuclear mechanism, the electron spin 
decay can be suppressed by app lying a _B-field or by po- 
larizing the nuclear system jl3l jlj| . 

In this Letter, we show that the spin T2 time, caused 
by spin-orbit interaction in GaAs QDs, is as large as the 
spin Ti time. We assume low temperatures, T <C Tiloq, 
where Tiujq is the size-quantization energy of the dot, and 
with no external noise in the applied magnetic field B. 
We, thus, argue that the lower bound Ti > 50/Lts es- 
tablished in Ref. 018, in fact, also a lower bound for T2. 
Furthermore, we show that the spin decay can be reduced 
by a special choice of direction of B, if there is Rashba 
spin-orbit interaction. 

The Hamiltonian describing the electron in a QD reads 
H = Hd + Hz + Hso + Uph, (1) 
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where p = —ih'V + (e/c)A(r) is the electron 2D kinetic 
momentum, U (r) is the lateral confining potential, with 
r = {x,y), and cr are the Pauli matrices. The axes x 
and y point along the main crystallographic directions 
in the (001) plane of GaAs. The spin-orbit coupling is 
given by the Hamiltonian (Q, where the term propor- 
tional to /3 originates from the bulk Dresselhaus spin- 
orbit coupling, which is due to the absence of inversion 
symmetry in the GaAs lattice; the term proportional to 
a represents the Rashba spin-orbit coupling, which can 
be present in quantum wells if the confining potential (in 
our case along the z-axis) is asymmetric. The magnetic 
field B = i3(sin6'cos(^,sin0sin(^,cos0) defines the spin 
quantization axis via the Zeeman term |0J. The phonon 
potential is given by 

U,,ir)^Y. ^l'^^^'T (^/3q.-»gSq,-)(&U,-+6q,-) (5) 
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where bl^ creates an acoustic phonon with wave vector 
q = (q||,(?2), branch index j, and dispersion ujqj; pc is 
the sample density [volume is set to unity in ©]. The 
factor F{ciz) equals unity for \qz\ ^ dr^ and vanishes 
for 1 ^ , where d is the characteristic size of the 
quantum well along the z-axis. We take into account 
both piezo-electric (/3qj ) and deformation potential (Sqj ) 
kinds of electron-phonon interaction. For /3qj and Sqj , we 
refer the reader to Ref. Next, we derive an effective 
Hamiltonian, which describes the spin dynamics, deco- 
herence and relaxation at low temperatures T <C Tiujq. 

The electron spin couples to phonons due to the spin- 
orbit interaction Q. For typical GaAs QDs the spin- 
orbit length \so = fi/iTi*(3 is much larger than the elec- 
tron orbit size A. In 2D, there is no linear in X/Xso 
contribution to the spin- pho non coupling at zero Zee- 
man sphtting [ifl lli llfl ll3 ■ Here, we consider a finite 
magnetic field, m*/3^ <C gpsB <C ?iwo, for which the lin- 
ear in A/ Ago contribution to the spin-phonon coupling 
dominates the spin decay. Using perturbation theory (or 
Schricffcr- Wolff transformation), we obtain the effective 
Hamiltonian Il9l 



^5Ms(B-)-5B(t)) - cr. 



(6) 
(7) 

(8) 



Here, is the electron orbital wave function, Ld is 
the dot Liouvillean, LdA = [Hd,A\. The vector ^ has 



HcS 

5^{t) = 2B X n{t), 
where ^{t) is defined as follows 



a simple form in the coordinate frame: x' = [x + y)/ \/2, 
V' = -{x-y)/V2, z' = z, namely, | = {y' /X-,x' /\+,0), 
where 1/A± = TO*(/3±a)/?i. Eq. {T)) contains one of our 
main results: in 1st order in spin- orbit interaction, there 
can be only transverse fluctuations of the effective mag- 
netic field, i.e. (5B(i) • B = 0. This statement holds true 
for spin coupling to any fluctuations, be it the noise of a 
gate voltage or coupling to particle-hole excitations in a 
Fermi sea. Next, we consider the decay of the electron 
spin, S = <t/2, governed by the Hamiltonian ©. 

The phonons which are emitted or absorbed by the 
electron leave the dot during a time Tc, d/s ^ Tc ^ X/s, 
where s is the sound velocity. The electron spin de- 
cays over a much longer time span in typical structures, 
and, therefore, undergoes many uncorrelated scattering 
events. In this regime, the dynamics and decay of the 
spin is governed by the Bloch equation 



(s) =LJx (s) -r(s) 



(9) 



where u; = ujl, with lu = gpBB/fi and I = B/_B. For 
a generic (5B(i}, we find that in the Born-Markov ap- 



proximation |2Cl |2] 

r = F'' + F'^, with 



the tensor Fi, can be written as 



F[j = 5ij {Spq - Iplq) J+ (w) - (S^p - lilp)J+{uj) 

-SijEkpqhlpqi^) + eipqlplgj{u}), (10) 

rf, = s,,ipiqj+{o) - hipj+{o), (11) 

where J^{w) = Ile[Jij{w) ± Jij{—w)] and I^{w) — 
Im \ Jij(vj) ± Jij{—vj)\ are given by the spectral function 



= ^J^ {SB,iO)6B,it)) e-^-'dt. (12) 
The inhomogeneous part in Eq. jnj is given by 



2Ti = ljJ,jiuj)-hJjjiuj)+eipqI+{uj) 

+e^qklklp [lpq{i^) - IpqiO)] , 



(13) 



where Sijk is the anti-symmetric tensor (with Ein- 
stein summation convention) and we have assumed that 
((5B(t)) = 0. Eq. @ describes spin decay in a number 
of problems, such as electron scatterin g off impurities in 
bulk systems, nuclear spin scattering |20j, etc. In our 
notation, the spin decay comes from the symmetric part 
of F, whereas the anti-symmetric part leads to a correc- 
tion to u) in Eq. ©. The tensor F*" describes spin decay 
due to processes of energy relaxation, such as, e.g., emis- 
sion/absorption of a phonon. Therefore, the Ti time is 
entirely determined by F*", see below. The tensor F"* can 
be non-zero only due to elastic scattering of spin, i.e. due 
to dephasing. F'' contributes to the decoherence time T2, 
and so does F*". In many cases, however, the latter con- 
tribution is negligible, and F'' entirely dominates the spin 
decoherence j22|. This is in strong contrast to what we 
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FIG. 1: Solid curve: the relaxation rate 1/Ti, given in 
Eq. II2H . as a function of the B-field for a GaAs QD with 
huJo = 1.1 meV, Ago = Ti/m*l3 — 1 fim, and a = 0. The di- 
rection of B is in plane, B = {Bx, By, 0). The dashed (dotted) 
curve shows the contribution of the transverse (longitudinal) 
phonons interacting with the electron via the piezo-electric 
(Pj^) mechanism. The dot-dashed curve shows the contribu- 
tion of the longitudinal phonons interacting via the deforma- 
tion potential (Ej) mechanism. The angular dependence of 
l/Ti is given by Eq. (1221 . the square of an ellipsoid. 



find here for an electron localized in a QD. To illustrate 
this, we first consider an example when dominates the 
spin decoherence and then return to our case. A textbook 
example is {SBi{0)5Bj{t)) = b'^Sij cxp{—\t\/Tc). Choos- 
ing I = (0, 0, 1), we obtain from Eqs. ifTUI) - ifT^ the non- 
zero elements: T^, = T^y = TlJ2 = jfPr,/ (1 + lu^t^), 
and V^^ = Tyy = 7^&^Tc, where 7„ = gtJ-B/'h. The longi- 
tudinal component {Sz) decays over the time Ti = = 
The transverse components decay over the time 
T2 = l/iTl^+TiJ. At w > 1/tc, the contribution of 
to T2 is negligible, and hence, T2 ^ Ti. The latter rela- 
tion has widely been quoted in the literature on quantum 
computation. In stark contrast to this example, we show 
now below that there are no intrinsic dephasing mecha- 
nisms for our case, which would justify this relation for 
the electron spin in GaAs QDs at T ^ hujo- 

We start with calculating the spin decay due to the 
mechanism Q. Here, Tfj is identically zero, due to the 
transverse nature of the fluctuating field (5B. This can 
be inferred from Eqs. (|ll|l and H12|l . noticing that each 
term in pi(l contains I ■ S'B — 0. In order to calculate 
T^j, we first find the main axes of the tensor Jij{w), see 
Eq. 1)12(1 . Jij{w) is diagonal in the frame {X, Y, Z), which 
is obtained from (a;', y' , z) by a rotation with Euler angles 
6*, and x- Here, the angles ip' = ip — it/ A and 9 give B 
in the frame (a;', y' , z), and x depends on the details of the 
dot confining potential U{r). It can be determined from 
(SBxSBYit)) = 0. Assuming C/(r) = U{r), we find [H 
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We now consider a harmonic confinement, U{r) = 
TO*WQr^/2, and evaluate fl{t) of Eq. ((SJ for the groimd 
state tp{r) — exp (— r^/2A^) /X^/n, where A~^ = 
h'^y/{m*u;o)'^ + (eB^/2c)2. Using the identity |2^ 
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we find an expression for ^l^' , which can be obtained from 
the r.h.s. of Eq. by replacing 



exp (iqii 



2iqy 



m*Ci;gA_ 



exp (-gj|AV4 



(16) 



Furthermore, an expression for is obtained in the 
same way, using the substitution (|16(l with a different 
prefactor, namely with (g^'/A-) {qx'/X+)- Finally, 
for the real part of the correlator Jij{w), we obtain 



ReJxx{w) 



sin^ § 




where — {e!^^l^ — l) , and Sj is the souird velocity 
for branch j. For GaAs, we use s\ ~ 4.7 x lO^cm/s and 
S2 = S3 ~ 3.37 X 10^ cm/s. Furthermore, = (5j,iE!o 
with So ~ 1 cV , and f-i\{) — 3\/27r/ii4K~^ sin^ iScosi?, 
(^2^) — \/27r/ii4K^^ sin2i9, /Sg^ = V27r/ii4K^-^(3 cos^ — 
1) sin I?, with ^,14 w 0.16 C/m^ and k ~ 13. The effective 
spin-orbit length A+ in Eq. H17|) is given by 



1^ 

A^ 
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4^2 



A^A2 
(18) 



ReJYy(w) is obtained from Eq. (|17|l by substituting 
A+ ^ A_, and Jzziw) = 0. ImJxxiw) and ImJyy(u') 
are irrelevant for our discussion, see further. From 
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Eq. (dni), we obtain 

^\z = Jxx{^) + Jyy^'^^ 

l^^iio). Since Fy/tj ^ tjA^/woAj. ^ 1, we can solve 
the secular equation iteratively and obtain [2^ 

1 

Then, the solution of Eq. Q reads {Sx{t)) — 
5_Le-*/^^ sin(cjt -t- </>), {SY[t)) = S'j.e"*/^^ cos(wi + 0), 
and {Sz{t)) ^ St + (S'| - S't) e"*/^! , with the ther- 
modynamic value of spin being St — l{l ■ T) Ti ~ 
-{l/2)ta,nh{hu;/2kBT), and the initial value (S(0)) = 
(5'_L sin^, S'_L cos(/>, 5^). For our special situation with 
purely transverse fluctuations (F'^ = 0), we obtain 

1 2 



Ti To 



Jxxi^) + JyYi'^)- 



(21) 
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with dimen- 

- (^/c)^ = 1, 
^^2 



Keeping in mind the setup of Ref. 0, we plot the re- 
laxation rate 1/Ti as a function of B for 9 = tt/2 and 
a ~ on Fig. n(sohd curve). We find that l/Tj has a 
plateau in a wide range of magnetic fields (cf. Ref.|3j), due 
to a crossover from the piezoelectric-transverse (dashed 
curve) to the deformation potential (dot-dashed curve) 
mechanism of electron-phonon interaction. For arbitrary 
(fi, e, and a, we have l/Ti = f /Ti[e = 0, a = 0) with 

/ = -L \{a^ + /?2)(i + cos^ 9) + 2al3s\xv^ 9s\xv2iA (22) 

The angular dependence of 1 /Ti for a = has been ob- 
tained before . Note that ^/J describes an ellipsoid in 
the frame (x', y' , z), i.e. f = x' + y' + 
sionless x',y',z obeying (x'/a)^ + {y' /h)"^ 
where a = 1 + a//?, 6=1 — a/ 13, and c = vt?+T?- 
Note that, if a = /3, then 6 = 0, i.e. 1/Ti vanishes 
if the B- field is along y'. The same is true along the 
x' axis for a = —(3. This special case (a = ±/3) has 
been discussed previously for extended electron states in 
a two-dimensional electron gas (2DEG) H^. Note that 
the Hamiltonian conserves the spin component <Jyi(x') 
for a = /3 (a = —(5) and B || y' {x'). This spin conser- 
vation results in Ti being infinite to all orders in the 
spin-orbit Hamiltonian At the same time, the deco- 
herence rate I/T2 reduces to the next order contribution 
of ©• However, as we show below, a single-phonon pro- 
cess is inefficient in inducing dephasing and, therefore, 
1 /T2 can be non-zero only in the next order in electron- 
phonon interaction. Next, we note that a long lived spin 
state also occurs in a different GaAs structure, namely, 
for a 2DEG grown in the (110) crystallographic direction. 
Then, the normal to the 2DEG plane component of spin 
is conserved [2^, provided there is no Rashba coupling. 

We discuss now other spin-orbit mechanisms. In 
Eq. we did not include the so-called A:'^-terms of 
the Dresselhaus spin-orbit coupling |2^, i.e. Hso oc 
f3(P {^(Jx{kxj ky} — cTy{ky, kx})- They are parametrically 
small {(P/X^ <C 1) in the 2D hmit, compared to the re- 
tained ones. However, their contribution to the spin de- 
cay can be important, if 5 ^ d^/A^ and <C dXso, 
since the orbital effect of the magnetic field contributes 
here in the first place. Still, for typical GaAs QDs, these 
mechanisms are negligible. 

Additional spin decay mechanisms arise from the di- 
rect spin-phonon interaction [loj . The strain field pro- 
duced by phonons couples to the electron spin via the 
spin-orbit interaction, resulting in the term AH' = 
{Vo / 4:)eijkO'i{uij , pk} , where pi is the bulk kinetic mo- 
mentum, Uij is the phonon strain tensor, and Vq — 
8 X lO^cm/s for GaAs. A similar mechanism occurs in 
a i?-field, due to g-factor fluctuations caused by lattice 
distortion. This yields AH" = gfiB UijaiBj, where 
g 10 for GaAs. The contribution of these mecha- 
nisms to the spin-flip rates in QDs has been estimated 
in Ref. ^3- Except for the a = ±/3 cases discussed 



above, the direct mechanisms are usually negligible in 
QDs. Here, we find that such spin-phonon couplings do 
not violate the equality T2 = 2Ti [23. For this, we note 

that Tfj = for a generic SBi = Y,q Mi{q){b^_^ + 5q) in 
2 



Eq. ®, if g |Af(q)|^ ^ at g ^ 0. Obviously, this con- 
dition is satisfied for the direct spin-phonon mechanisms, 
since the strain tensor Uij vanishes at q = 0. The same 
follows for the Hamiltonian |^ with the phonon potential 
Q and an arbitrary Hso ! the physical explanation being 
that the potential of long-wave phonons is constant over 
the dot size and, thus, commutes with Hso- Finally, 
we note that, at temperatures T ~ hujo, there can be 
dephasing mechanisms |28l | , which can result in T2 <C Ti . 

In conclusion, we have shown that the dccoherence 
time T2 of an electron spin in a GaAs QD is as large 
as the relaxation time Ti for the spin decay based on 
spin-orbit mechanisms. We acknowledge support from 
the Swiss NSF, NCCR Basel, EU RTN 'Spintronics', US 
DARPA, ARO, and ONR. 
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